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Abstract — In  an  irredundant  sum-of-products  expression  (ISOP),  each  product  is  a  prime  impiicant  (Pi)  and  no  product  can  be  deieted 
without  changing  the  function.  Among  the  ISOPs  for  some  function  /,  a  worst  ISOP  (WSOP)  is  an  ISOP  with  the  iargest  number  of  Pis 
and  a  minimum  iSOP  (MSOP)  is  one  with  the  smaiiest  number.  We  show  a  class  of  functions  for  which  the  Minato-Morreaie  iSOP 
aigorithm  produces  WSOPs.  Since  the  ratio  of  the  size  of  the  WSOP  to  the  size  of  the  MSOP  is  arbitrariiy  iarge  when  n,  the  number  of 
variabies,  is  unbounded,  the  Minato-Morreaie  aigorithm  can  produce  resuits  that  are  very  far  from  minimum.  We  present  a  ciass  of 
muitipie-output  functions  whose  WSOP  size  is  aiso  much  iarger  than  its  MSOP  size.  For  a  set  of  benchmark  functions,  we  show  the 
distribution  of  iSOPs  to  the  number  of  Pis.  Among  this  set  are  functions  where  the  MSOPs  have  aimost  as  many  Pis  as  do  the 
WSOPs.  These  functions  are  known  to  be  easy  to  minimize.  Aiso,  there  are  benchmark  functions  where  the  fraction  of  iSOPs  that  are 
MSOPs  is  smaii  and  MSOPs  have  many  fewer  Pis  than  the  WSOPs.  Such  functions  are  known  to  be  hard  to  minimize.  For  one  ciass  of 
functions,  we  show  that  the  fraction  of  iSOPs  that  are  MSOPs  approaches  0  as  n  approaches  infinity,  suggesting  that  such  functions 
are  hard  to  minimize. 

Index  Terms — Logic  minimization,  compiete  sum-of-products  expressions,  irredundant  sum-of-products,  muitipie-output  functions, 
heuristic  minimization,  prime  impiicants,  symmetric  functions,  minimum  sum-of-products  expressions,  worst  sum-of-products 
expressions,  graph  enumeration,  minimaiiy  strongiy  connected  digraphs. 

-  -f  - 


1  Introduction 

wo-LEVEL  logic  minimization  is  a  basic  problem  in  logic 
synthesis.  Although  algorithms  exist  that  obtain  the 
exact  minimum  sum-of-products  expressions  (MSOP)  for  a 
large  set  of  functions  [7],  [8],  the  majority  of  practical 
systems  use  heuristic  logic  minimization  algorithms.  These 
produce  irredundant  sum-of-products  expressions  (ISOPs) 
that  are  not  necessarily  minimum.  For  example,  PRESTO 
[4],  [33],  MINI  [15],  ESPRESSO  [3],  and  others  [10],  [24] 
produce  nonmmimum  ISOPs. 

An  ISOP  is  the  OR  of  prime  impiicants  (Pis)  such  that 
deleting  any  PI  changes  the  function.  For  example,  two 
expressions  XiX2  V  X2X3  V  xiX's  and  X1X2  V  XiX^  V  X1X2  V 
X1X3  are  both  ISOPs  for  the  same  function  (see  Fig.  la  and 
Fig.  lb).  The  first  is  an  MSOP  and  the  second  is  a  worst 
ISOP  (WSOP),  an  ISOP  with  the  largest  number  of  Pis.  Most 
practical  logic  synthesis  algorithms  generate  ISOPs  at  some 
point,  so  an  understanding  of  ISOPs  is  crucial. 

In  this  paper,  we  show  classes  of  fcmctions,  where  the 
ratio  of  the  WSOP  size  (number  of  Pis)  to  the  MSOP  size  is 
arbitrarily  large  when  the  number  of  variables  is  un¬ 
bounded.  We  show  that  the  Minato-Morreaie  algorithm 
[19],  [20]  produces  WSOPs  for  this  class.  We  also  show  an 
M-variable  multiple-output  function  whose  MSOP  size  is  at 
most  2n  and  whose  WSOP  size  is  at  least  2". 
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We  also  show  an  algorithm  that  produces  all  ISOPs  for  a 
given  function.  When  applied  to  benchmark  functions,  we 
notice  a  correlation  between  the  degree  of  difficulty  in 
determining  an  MSOP  and  the  distribution  of  ISOPs  to  the 
number  of  Pis.  For  example,  easily  minimized  functions 
tend  to  have  a  larger  proportion  of  ISOPs  that  are  MSOPs. 

We  also  show  an  analysis  of  a  class  of  functions  with 
respect  to  the  number  of  MSOPs  and  WSOPs.  For  this 
class,  the  number  of  ISOPs  that  are  MSOPs  and  WSOPs  is 
n[n  —  1) !  and  respectively.  Since  there  are  many  more 
WSOPs  than  MSOPs,  it  suggests  that  such  functions  are 
difficult  to  minimize. 

This  paper  is  organized  as  follows:  Section  2  shows 
definitions  and  basic  properties.  Section  3  considers  the 
MSOPs  and  WSOPs  of  specific  functions.  The  ratio  of  the 
number  of  Pis  in  a  WSOP  to  that  of  the  MSOP  is  arbitrarily 
large  when  n  is  unbounded  for  these  functions.  Section  4 
presents  a  class  of  multiple-output  functions  for  which 
there  is  a  large  disparity  between  the  number  of  Pis  in  an 
MSOP  and  in  a  WSOP.  Section  5  focuses  on  the  distribution 
of  the  number  of  ISOPs  to  the  number  of  Pis  required.  It 
presents  a  method  to  derive  all  ISOPs  of  a  given  function. 
Section  6  shows  experimental  results  obtained  from  this 
method.  Section  7  shows  how  this  distribution  can  be 
obtained  analytically  for  specific  functions.  Section  8 
concludes  the  paper. 

2  Definitions  and  Fundamental  Results 

In  the  discussions  to  follow,  we  will  often  use  symmetric 
functions. 
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Fig.  1 .  Karnaugh  maps  for  5T(3, 1).  (a)  MSOP.  (b)  WSOP. 


Lemma  2.1. 


t{WSOP  :  5[o,i,3,4,6,7})  =  70. 

Proof.  See  the  Appendix.  □ 

3  WSOPs  AND  MSOPs  FOR  Specific  Functions 

Definition  3.1.  Let  ST{n,  k)  he  a  symmetric  function  of 
n-variables  xi,X2,  ■  ■  ■  ,Xn  such  that 


Definition  2.1.  S\,  a  (totally)  symmetric  function,  is  1  ifm  of 

its  n  variables  are  1,  where  m  e  A  C  {1, 2, . . . ,  n}  and  is  0 
otherwise. 

Example  2.1.  The  AND  and  OR  of  n  variables  are  symmetric 
functions,  represented  by  and  S'jfj  2  respectively. 

Definition  2.2.  x  and  x  are  literals  of  a  variable  x.  The  AND  of 
literals  is  a  product.  The  OR  of  products  is  a  sum-of- 
products  expression  (SOP). 

Definition  2.3.  A  prime  implicant  (PI)  of  a  function  f  is  a 
product  that  implies  f  such  that  the  deletion  of  any  literal  from 
the  product  results  in  a  new  product  that  does  not  imply  f. 

Definition  2.4.  A  complete  sum-of-products  expression 
(CSOP)  [2],  [22]  of  a  function  /  is  the  SOP  of  all  Pis  of  f. 

Definition  2.5.  An  irredundant  sum-of-products  expression 
(ISOP)  is  an  SOP  where  each  product  is  a  PI  and  no  PI  can  be 
deleted  without  changing  the  function  represented  by  the 
expression. 

Definition  2.6.  Among  the  ISOPsfor  f,  the  one  with  the  largest 
number  of  Pis  is  a  worst  ISOP  (WSOP)  and  the  one  with  the 
smallest  number  of  Pis  is  a  minimum  SOP  (MSOP). 

Defintion  2.7.  The  size  of  an  SOP  is  the  number  of  Pis  in  the 
SOP.  The  size  of  a  CSOP,  WSOP,  and  MSOP  of  function  f  is 
denoted  as  t{CSOP  :  /),  t{WSOP  :  /),  and  t{MSOP  :  f), 
respectively. 

The  following  is  well  known. 

Theorem  2.1  [13],  [22].  For  any  switching  function  of 
n  variables,  t{MSOP  :  f)  < 

This  upper  bound  is  firm.  For  example,  the  exclusive  OR 
function,  /exor  =  aii  ©  2:2  ©  . . .  ©  x„,  has  2"“^  minterms,  all 
of  which  are  Pis,  and,  so,  t{MSOP  :  /exor)  =  2"“^ 

Further,  t{WSOP  :  /exor)  =  2"“^  there  being  only  one 
ISOP.  It  is  tempting  to  believe  that  t{WSOP  :  /)  <  2"“^,  for 
any  f.  Indeed,  Meo  [18]  conjectured  this  in  the  mid  1960s. 
However,  a  counterexample  was  published  in  Russian  in 
1962  by  Yablonski  [37]  (which  was  reported  in  English  by 
Kautz  [16]  in  1966).  Specifically,  Yablonski  showed: 

Theorem  2.2  [37].  There  exists  a  switching  function  on 
n  variables  where  t{WSOP  :  f)  >  2"“^. 

by  showing  an  ISOP  for  S'jy  13457}  with  70  Pis.  This  is  six 
more  than  the  upper  bound  of  2^“^  =  64  as  conjectured  by 
Meo.  As  we  have  not  seen  a  proof  of  this  in  English,  we 
include  one  here.  We  extend  Yablonski's  result  by  showing 
that  his  ISOP  is  a  WSOP. 


ST(n,  k) 


_  j  1  k  <  Xi  <n  —  k 


1=1 

0  otherwise, 


where  Xi  is  the  number  of  variables  that  are  1  and 
n  >  2k. 

Example  3.1.  ST(n,  0)  =  1.  ST(n,  |),  for  even  n,  is  the  OR 
of  all  minterms  with  exactly  half  of  the  variables 
complemented. 

Lemma  3.1.  ST{n,  k)  can  be  represented  as 


Example  3.2. 

5T(n,l)  =  5^1,2,,.., „}  5^0,1,...,.-!} 

=  (ail  V  2:2  V  •  •  •  V  a;„)(:ri  V  *2  V  •  •  •  V  :r„) 


and 


=  (a;ia;2  V  a;ia;3  V  •  •  •  V  2;„_ia;„) 
(:Ei:r2  V  :ri:r3  V  •  •  •  V  :r„_i:r„). 


We  are  interested  in  the  sizes  of  the  CSOP,  an  MSOP,  and  a 
WSOP  for  ST(n,k).  Voight  and  Wegener  [36]  consider  the 
CSOP  and  MSOP  sizes  for  general  symmetric  functions, 
stating  expressions,  and  outlining  a  proof.  Our  next  result 
gives  the  CSOP,  MSOP,  and  WSOP  sizes  for  ST(n,k) 
functions.  A  complete  proof  is  given  in  the  Appendix. 

Theorem  3.1. 

1.  riCSOP  :  STin,  k))  = 

2.  t{M  SOP -.STin, k))={f). 

3.  t{WSOP  :  ST(n,  fc))  >  2  (")  -  iff) . 

Proof.  See  the  Appendix.  □ 

In  the  proof  of  Lemma  2.1,  we  showed  that 
riWSOP  :  ST(7, 3))  =  56,  which  is  six  more  than  the  lower 
bound  given  in  3  of  Theorem  3.1.  Therefore,  for  n  =  7  and 
fc  =  3,  the  lower  bound  is  not  tight.  However,  for  fc  =  1,  the 
lower  bound  is  exact.  That  is. 

Theorem  3.2. 


t{WSOP  :  ST{n,  1))  =  2n  -  2. 

Proof.  See  the  proof  of  Theorem  7.3  in  the  Appendix.  □ 
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A  special  case  of  these  theorems  occurs  when  n  =  3  and 
fc=  1. 

Example  3.3.  5'T(3, 1)  =  {xi  \/  X3){xi  V  2:2  V  2:3)  has  the 

following  properties: 

1.  t{CSOP  :  ST{3,1))  =  6. 

2.  t{MSOP  :  ST{3, 1))  =  3. 

3.  t{WSOP  :  ST{3,1))  =  4. 

Fig.  la  and  Fig.  lb  show  the  MSOP  and  WSOP  of 
ST{3, 1),  respectively.  Interestingly,  the  ISOP  generator 
of  Minato  [19],  which  is  based  on  Morreale's  [20] 
algorithm  produces  a  WSOP  for  ST{3, 1)  instead  of  an 
MSOP.  This  will  be  discussed  in  more  detail  later. 

Definition  3.2.  The  redundancy  ratio  of  a  function  f  is 

_  rjWSOP  :  /) 
t{MSOP  :  f)  ■ 

The  normalized  redundancy  ratio  of  an  n-variable 
function  f  is 

c^if)  =  VKf), 

where  t{WSOP  :  f)  and  t{MSOP  :  f)  are  the  sizes  of 
WSOPs  and  MSOPs.  If  this  ratio  is  small,  any  logic 
minimization  algorithm  will  do  well  since,  even  if  a  WSOP 
is  generated,  it  is  not  much  worse  than  an  MSOP.  On  the 
other  hand,  a  large  ratio  suggests  that  care  should  be  exercised. 
The  normalized  redundancy  ratio  is  normalized  with  respect  to 
the  number  of  variables.  It  is  a  convenience;  it  allows  one  to 
compare  the  redundancy  ratio  of  two  functions  with  a  different 
number  of  variables. 

From  the  expressions  for  t{MSOP  :  ST{n,k))  and 
t{WSOP  :  ST{n,  k))  given  in  Theorem  3.1,  we  can  state: 

Theorem  3.3. 


TABLE  1 

p  and  a  for  ST{n,  1)  versus  n 


n 

p 

(7 

2 

1.0000 

1.0000 

3 

1.3333 

1.1006 

4 

1.5000 

1.1066 

5 

1.6000 

1.0986 

6 

1.6667 

1.0889 

7 

1.7143 

1.0800 

8 

1.7500 

1.0725 

^  OC' 

^  2 

^  1 

exist  functions  where  p  is  larger  than  2.  Indeed,  we  show  a 
class  of  functions  in  which  p  increases  without  bound  as  n 
increases.  This  has  important  consequences  for  heuristics 
that  produce  ISOPs.  For  such  heuristics  there  is  the  prospect 
of  generating  an  ISOP  whose  size  is  much  larger  than  the 
minimum.  We  consider  this  topic  now. 

Definition  3.3.  Let  ST(m,ky  be  the  n  =  m  ■  r-variable 
function 

STirn.,  k^  (2:1 , 2:2 1  i  Xjyit) 

r 

STirn^  k)  1)+1 ) 

where  is  the  AND  (product)  of  r  functions. 

Theorem  3.5.  5T(m,  kf  has  the  following  properties: 

1.  r{CSOP  :  ST{m,  kf)  = 

2.  t{MSOP  :  ST(m,  kf)  =  (7)’'. 

3.  t{WSOP  :  ST{m,  kf  )  >  [2(™)  - 

Proof.  See  the  Appendix.  □ 


p{ST{n,  k))  >  2 


a{ST{n,  k))  > 


From  the  expressions  for  t{MSOP  :  ST{n,  1))  and 
t{WSOP  :  ST{n,  1))  given  in  Theorems  3.1  and  3.2,  respec¬ 
tively,  we  can  state: 

Theorem  3.4. 


For  fc  =  1,  we  have: 

Theorem  3.6. 

t{WSOP  :  ST{m,  I)")  =  2'’(m  -  1)’'. 

Example  3.4.  For  m  =  3  and  k  =  1,  ST{3,  l)"^  has  Pis, 
t{MSOP  :  ST{3,  !)'■)  =  3’’,  and 

t{WSOP  :  ST{3, 1)")  =  AP 

We  have: 


p{ST{n,  1))  = 
a{ST{n,  1))  = 

Table  1  shows  the  values  of  p  and  a  for  ST{n,  1),  where  p{ST{m,  1)'^)  =  2’’ 

2  <  n  <  8.  It  can  be  seen  that  a  takes  its  maximum  value 
when  n  =  4.  That  is,  as  n  increases  above  2,  ^  first  increases,  3.5.  por  m  =  4  and  fc  =  1,  we  have 

peaking  at  4,  and  then  it  continually  decreases. 

From  Theorem  3.4,  p  is  monotone  increasing  with  an  p(ST(A,  l)”^)  =  (l.S)*^. 

upper  limit  of  2.  Thus,  for  ST{n,  1)  functions,  the  number  of 

Pis  in  a  WSOP  is  never  more  than  two  times  the  number  of  From  this,  it  can  be  seen  that  p  becomes  arbitrarily  large  as  r 
Pis  in  an  MSOP.  An  important  question  is  whether  there  approaches  infinity.  In  this  example,  there  are  n  =  4r 
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variables.  This  represents  a  class  of  functions  for  which  p 
grows  without  bound  as  the  number  of  variables  grows. 

4  Extension  to  Multiple-Output  Functions 

In  the  case  of  multiple-output  functions,  minimization  of 
AND-OR  two-level  networks  or  programmable  logic 
arrays  (PLAs)  can  be  done  using  characteristic  functions 
[26],  [27],  [29]. 

Definition  4.1.  for  an  n-variable  function  with  m  output  values, 

fj{Xl,X2,...,Xn)  (j=  0,  1), 

form  an  {n  +  l)-variable  two-valued  single  output  func¬ 
tion  F{xi,X2,  ■  ■  ■  ,Xn,  X„+i),  where  Xi  is  a  binary  valued 
variable,  for  I  <  i  <  n  and  A„+i  takes  m  values  such 
that  F{xi,X2,...,Xn,j)  =  I  iff  fj{Xl,X2,  ■  ■  ■  ,x„)  =  I 
(j  =  0, 1, . . . ,  m  —  1).  Then,  F  represents  all  and  only  the 
permitted  combinations  of  inputs  and  nonzero  output  values  of 
f.  F  is  called  the  characteristic  function  (for  nonzero 
outputs). 

The  significance  of  the  characteristic  function  is  seen  in 
Theorem  4.1  below. 

Definition  4.2.  is  a  literal,  where  X  takes  a  value  in 
{0, 1, ...  ,p  —  1}  and  S'  C  {0, 1, ...  ,p  —  1}  such  that  =  1 
if  X  =  a  e  S  and  X^  =  0,  otherwise.  A  logical  product  of 
literals  that  contains  at  most  one  literal  for  each  variable  is  a 
product  term.  Products  combined  with  OR  operators  form  a 
sum-of-products  expression  (SOP).  A  Prime  implicant  (PI), 
irredundant  sum-of-products  expression  (ISOP),  worst  ISOP 
(WSOP),  and  minimum  SOP  (MSOP)  are  defined  in  a 
manner  similar  to  the  two-valued  case. 

Theorem  4.1  [15],  [27],  [29].  The  number  of  AND  gates  in  the 
minimum  AND-OR  two-level  network  for  the  function 
(/o;  /i;  •  •  ■  j  fm-i)  Is  cqual  to  the  number  of  Pis  in  the  MSOP 
for  the  characteristic  function  F. 

Definition  4.3.  An  n-bit  decoder  has  n  inputs  xi,X2, . .  ■  ,x„, 
and  2"  outputs  /o,  /i, . . . ,  /2"-i,  where  fi  =  0iff  the  binary 
number  representation  of  xiX2,  ■ . .  ,x„  is  i. 

Example  4.1.  The  4-bit  decoder  has  16  outputs,  as  follows: 

fo  =  XiV  X2V  X3V  Xi, 
fi  =  XiV  X2V  X3V  Xi, 

/i5  =  Xi  V  :f2  V  :r3  V  Xi. 

Definition  4.4.  DEC_n  is  the  characteristic  function  of  an  n-bit 
decoder. 

Example  4.2.  DEC_4  is  shown  in  positional  cube  notation  in 
the  upper  table  of  Fig.  2.  That  is,  each  entry  in  this  table  is 
a  prime  implicant  of  DEC_4,  where  Xi  appears  as  Xi,  Xi, 
or  don't  care  (absent)  if  the  corresponding  entry  is  10,  01, 
or  11,  respectively.  For  X5,  the  entry  0111111111111111  is 
the  literal  etc.  Therefore,  the 

first  entry,  10  10  10  10  0111111111111111,  corresponds  to 


ISOPs  with  16  Pis 


Xl 

X2 

X3 

Xi 

A5 

01 

01 

01 

01 

01  •••  15 

10 

10 

10 

10 

0111111111111111 

10 

10 

10 

01 

1011111111111111 

10 

10 

01 

10 

1101111111111111 

10 

10 

01 

01 

1110111111111111 

10 

01 

10 

10 

1111011111111111 

10 

01 

10 

01 

1111101111111111 

10 

01 

01 

10 

1111110111111111 

10 

01 

01 

01 

1111111011111111 

01 

10 

10 

10 

1111111101111111 

01 

10 

10 

01 

1111111110111111 

01 

10 

01 

10 

1111111111011111 

01 

10 

01 

01 

1111111111101111 

01 

01 

10 

10 

1111111111110111 

01 

01 

10 

01 

1111111111111011 

01 

01 

01 

10 

1111111111111101 

01 

01 

01 

01 

1111111111111110 

ISOPs  with 

8  Pis 

Xl 

X2 

X3 

Xi 

A5 

01 

01 

01 

01 

01  ■■ ■  15 

10 

11 

11 

11 

0000000011111111 

01 

11 

11 

11 

1111111100000000 

11 

10 

11 

11 

0000111100001111 

11 

01 

11 

11 

1111000011110000 

11 

11 

10 

11 

0011001100110011 

11 

11 

01 

11 

1100110011001100 

11 

11 

11 

10 

0101010101010101 

11 

11 

11 

01 

1010101010101010 

Fig.  2.  Positional  cubes  for  two  ISOPs  of  DEC_4. 

,,  .  .  ,  -  -  -  -  2.3,4, 5,6.7,8,9.10.11. 12, 13, 14,15} 

the  prime  implicant  XiX2X3XiX^ 

Collectively,  the  16  entries  in  the  upper  table  of  Fig.  2 
represent  an  ISOP  of  DEC_4  with  16  Pis.  An  ISOP  for 
DEC_4  with  only  eight  Pis  exists,  as  shown  in  the  lower 
table  of  Fig.  2. 

The  observations  of  Example  4.2  can  be  generalized  as 
follows: 

Theorem  4.2.  The  function  DEC_n  has  a  WSOP  that  requires  at 
least  2"  Pis  and  an  MSOP  that  requires  at  most  2n  Pis. 

The  above  theorem  proves  the  existence  of  an  n-variable 
2"-output  variable  function,  where  the  sizes  of  the  MSOP 
and  the  WSOP  are  at  most  2n  and  at  least  2",  respectively. 
The  upper  ISOP  for  DEC_4  shown  in  Fig.  2  is  not  a  WSOP 
since  an  ISOP  with  20  Pis  has  been  found  for  DEC_4. 

5  Derivation  of  All  ISOPs 

Very  little  is  known  about  the  distribution  of  the  sizes  for 
ISOPs.  For  example,  even  for  single-output  functions,  we 
know  of  no  study  that  shows  how  many  ISOPs  exist  with 
various  number  of  product  terms. 

Although  various  methods  to  generate  all  the  ISOPs  for  a 
logic  function  are  known  [22],  [12],  [21],  [6],  [35],  [25],  no 
experimental  results  have  been  reported.  Experiments  are 
computationally  intensive  even  for  functions  with  a  small 
number  of  variables.  However,  we  can  obtain  the  statistical 
properties  of  ISOPs  for  some  interesting  functions. 

Before  showing  the  complete  algorithm,  consider  the 
following: 
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Fig.  3.  f(xi,X2,X3)  =  ST{3,1). 


Example  5.1. 

/  =  ST{3, 1)  =  {xi  V  012  V  X3){xi  V  ^2  V  xs) 

has  six  minterms  (Fig.  3)  and  six  Pis 

[Pi  :  ^12:3,  P2  ■  X1X2,  P3  ■■  X2X3,  Pi  :  2:1^3,  P5  :  2:1^2,  Pe  :  ^22:3]- 

Fig.  4  is  the  covering  table  for  ST{3, 1).  It  shows  the 
following  relations: 

To  cover  mi,pi  V  pe  is  necessary. 

To  cover  m^^pi  V  p2  is  necessary. 

To  cover  m2,p2  V  p3  is  necessary. 

To  cover  mQ,p3  V  pi  is  necessary. 

To  cover  mi,pi  V  ps  is  necessary. 

To  cover  m^^p^  V  pe  is  necessary. 

To  satisfy  all  the  conditions  at  the  same  time,  we  have 
P{f)  =  1,  where 

P{f  )  =  (Pl  V  P6)(pi  V  P2){P2  V  P3)(P3  V  P4)(P4  V  P5)(p5  V  pg)- 

P(/)  is  called  the  Petrick  function  [22].  By  expanding 
P{f)  into  SOPs,  we  have 

P{f)  =  (pi  V  P2P6)(P3  V  P2P4)(P5  V  PiPe) 

=  PlPiPh  V  P2P3P5P6  V  P1P2P4P5  V  P2P4P5P6 

V  P1P3P4P6  V  P2P3P4P6  V  P1P2P4P6  V  P2P4P6- 

Note  that  each  product  with  an  underline  is  covered  by 
another  product  having  fewer  literals.  Such  products  are 
redundant.  Deleting  these  products,  we  have 

P{f)  =  P1P3P5  V  P2P3P5P6  V  P1P2P4P5  V  piP3P4Pe 
V  P2P3P4P6  V  P2P4P6- 

P(/)  consists  of  all  the  Pis  of  the  Petrick  function  [23]  and 
each  PI  of  P(/)  corresponds  to  an  ISOP  for  /. 
Furthermore,  each  literal  p;  in  the  PI  of  P(/)  corresponds 
to  a  PI  for  /.  For  example,  pipsps  corresponds  to  the  ISOP 
xiXz  V  a;22;3  V  a:ia:2-  Note  that  there  are  six  ISOPs;  two 
have  three  Pis,  while  four  have  four  Pis.  Thus,  5'T(3, 1) 


minterms 


m\ 

in-i 

m2 

nifi 

777,5 

Pl 

1 

1 

P2 

1 

1 

Pi 

1 

1 

Pis  Pi 

1 

1 

PS. 

1 

1 

PC, 

1 

1 

has  two  MSOPs  with  three  Pis  and  four  WSOPs  with 
four  Pis. 


In  this  way,  all  the  ISOPs  are  obtained.  For  general 
functions,  the  number  of  minterms  and  Pis  are  very  large. 
Thus,  we  use  an  ROBDD  (reduced  ordered  binary  decision 
diagram)  to  represent  the  function  and  a  Prime_TDD 
(Ternary  decision  diagram)  [31]  to  represent  the  set  of  all 
the  Pis.  In  the  Prime_TDD  for  /,  each  path  from  the  root 
node  to  the  constant  1  node  corresponds  to  a  PI  for  /.  We 
also  use  an  ROBDD  to  represent  the  Petrick  function.  While 
there  are  many  ways  to  generate  all  the  ISOPs  of  a  given 
function  /,  we  use  the  following  algorithm: 

Algorithm  5.1  (Generation  of  all  ISOPs  for  a  function  /). 

1.  Generate  all  the  Pis  for  /  by  using  the  Prime_TDD 
(the  ternary  decision  diagram  representing  Pis)  of  /. 

2.  From  the  set  of  Pis  and  the  set  of  minterms  for  /, 
generate  the  Petrick  function  P{f)  (which  represents 
the  covering  table  [22]). 

3.  Generate  the  Prime_TDD  (which  represents  all  the 
Pis)  of  P{f). 

4.  Generate  the  1-paths  of  the  Prime_TDD  and,  for  each 
1-path,  generate  the  corresponding  ISOP. 

In  the  Prime_TDD  in  Step  4,  each  path  from  the  root  node  to 
the  constant  1  corresponds  to  a  PI  for  P{f)  and  to  an  ISOP 
for  /.  Each  1  edge  has  weight  1  and  each  0  edge  has 
weight  0.  The  total  sum  of  weights  from  the  root  node  to  the 
constant  1  nodes  is  the  number  of  Pis  in  the  ISOP.  Note  that 
the  shortest  path  corresponds  to  an  MSOP  and  the  longest 
path  corresponds  to  a  WSOP. 

6  Experimental  Results 

6.1  ST{n,k)  Functions 

Using  Algorithm  5.1,  we  compare  the  number  of  Pis  in 
ST{n,  ky  for  different  n,  k,  and  r.  Table  2  shows  the  number 
of  Pis  in  the  MSOP  and  the  WSOP  of  ST(n,  ky,  as  well  as 
the  total  number  of  Pis.  Shown  also  are  the  results  of  the 
Minato-Morreale  algorithm. 

The  9SYM  (or  SYM9)  [11],  [15]  function  shown  in  [3, 
p.  165]  is  identical  to  ST(9,3).  It  has  1,680  Pis, 
t{WSOP  :  ST{9,  3))  >  148,  and  t{MSOP  :  ST{9,  3))  =  84. 
POP  [9],  a  PRESTO-type  [4],  [33]  logic  minimization 
algorithm,  produced  an  ISOP  with  148  products.  CAMP 
[1]  produced  an  ISOP  with  130  Pis,  while  MINI  [15]  did 
well,  producing  85  Pis. 

Table  3  shows  the  distribution  of  ISOPs  to  the  number  of 
Pis  in  an  ISOP  for  ST{n,  1)  for  3  <  n  <  7.  This  data  was 
obtained  by  Algorithm  5.1.  It  can  be  seen  that  the  set  of 
MSOPs  is  small  compared  to  the  set  of  all  ISOPs. 

6.2  Other  Functions 

We  also  applied  Algorithm  5.1  to  compare  the  number  of 
Pis  for  multiple-output  functions.  Table  4  shows  the 
distribution  of  the  number  of  Pis  in  ISOPs  for  various 
arithmetic  functions  [32]. 

INCn  is  an  n-input  n  +  I  output  function  such  that  the 
value  of  the  output  is  x  +  I,  where  x  is  the  value  of  the 


Fig.  4.  Covering  tabie  of  5T(3, 1). 
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TABLE  2 

Number  of  Pis  and  Redundancy  Ratio  for  Various  Functions 


n 

MSOP 

WSOP 

ISOP 

MM 

PI 

P 

a 

sr(3,i) 

3 

3 

4 

4 

6 

1.33333 

1.1006 

5r(3,i)^ 

6 

9 

16 

16 

36 

1.77778 

1.1006 

ST(3,lf 

9 

27 

64 

64 

216 

2.37037 

1.1006 

ST(3,lf 

12 

81 

256 

256 

1296 

3.16049 

1.1006 

ST(3, 1)5 

15 

243 

1024 

1024 

7776 

4.21399 

1.1006 

sr(3,i)‘^ 

18 

729 

4096 

4096 

46656 

5.61866 

1.1006 

57(4,1) 

4 

4 

6 

6 

12 

1.50000 

1.1067 

57(4,1)2 

8 

16 

36 

36 

144 

2.25000 

1.1067 

57(4,1)3 

12 

64 

216 

216 

1728 

3.37500 

1.1067 

57(4,1)3 

16 

256 

1296 

1296 

20736 

5.06250 

1.1067 

57(4, 1)5 

20 

1024 

7776 

7776 

248832 

7.59375 

1.1067 

57(5,1) 

5 

5 

8 

8 

20 

1.60000 

1.0986 

57(5,1)2 

10 

25 

64 

64 

400 

2.56000 

1.0986 

57(5,1)3 

15 

125 

512 

512 

8000 

4.09600 

1.0986 

57(5,1)3 

20 

625 

4096 

4096 

160000 

6.55360 

1.0986 

57(5,2) 

5 

10 

14 

14 

30 

1.40000 

1.0696 

57(6,1) 

6 

6 

10 

10 

30 

1.66667 

1.0889 

57(6,2) 

6 

15 

24 

24 

90 

1.60000 

1.0815 

57(7, 1) 

7 

7 

12 

12 

42 

1.71429 

1.0800 

57(7,2) 

7 

21 

t 

36 

210 

t 

t 

57(7,3) 

7 

35 

t 

50 

420 

t 

t 

57(8,1) 

8 

8 

14 

14 

56 

1.75000 

1.0725 

57(8,2) 

8 

28 

t 

50 

420 

t 

t 

57(9,1) 

9 

9 

16 

16 

72 

1.77778 

1.0660 

57(9,2) 

9 

36 

t 

66 

756 

t 

t 

57(10.1) 

10 

10 

18 

18 

90 

1.80000 

1.0605 

n:  nuiiiber  of  input  variables. 

MSOP:  number  of  Pis  in  MSOP. 

WSOP:  number  of  Pis  in  WSOP. 

ISOP  MM:  immber  of  Pis  in  an  ISOP  generated  by  Minato-Morreale’s 
algoritlnii. 

PI:  nmiiber  of  prime  implicants. 

p:  redundancy  ratio:  t{WSOP  :  f)/T{MSOP  :  /) 

a:  normalized  redundancy  ratio: 

f:  memory  overflow  precluded  generation  of  these  values 


input;  WGT5  is  the  same  as  RD53,  a  5-input  3-output 
function,  where  the  output  is  a  binary  number  whose  value 
is  the  number  of  Is  on  the  inputs;  ROT6  computes  the 
square  root  of  a  6-bit  integer;  LOG5  computes  the  logarithm 
of  the  5-bit  integer;  ADR3  is  a  3-bit  adder;  and  SQR5 
computes  the  square  of  the  5-bit  input. 

Note  that  all  the  ISOPs  for  INC6  have  the  same  number 
of  Pis.  This  means  any  logic  minimizer  obtains  an  exact 
minimum  solution.  This  is  also  is  true  for  WGT5.  For  ADR3, 
most  of  the  ISOPs  have  31  Pis  or  33  Pis.  This  is  consistent 
with  the  observation  that  the  logic  minimization  of  ADR3  is 
relatively  easy.  For  SQR5,  the  distribution  is  very  wide.  The 
MSOPs  have  27  Pis,  while  WSOPs  have  37  Pis.  This  is 
consistent  with  the  observation  that  the  minimization  of 
SQR5  is  more  difficult.  Note  that  SQR5  is  a  10  output  binary 
function.  The  data  shown  is  for  all  outputs. 

Although  we  could  not  obtain  the  distribution  for  SQR6 
due  to  the  memory  overflow,  we  conjecture  that  the 


distribution  of  number  of  Pis  for  SQR6  is  also  wide.  We 
also  developed  WIRR,  a  heuristic  algorithm  to  obtain  ISOPs 
with  many  products.  For  SQR5,  SQR6,  and  9SYM,  the 
numbers  of  Pis  in  the  solutions  are  shown  in  Table  5. 

7  Distribution  of  ISOPs — An  Analytic 
Approach 

The  distribution  of  ISOPs  to  the  number  of  Pis  is  a  way  to 
represent  the  search  space  a  heuristic  algorithm  must 
traverse  in  a  minimization  of  an  expression.  For  the  case 
of  ST{n,  1)  functions,  we  can  show  a  part  of  this  distribu¬ 
tion;  a  graph  representation  of  the  set  of  Pis  allows  this. 
Definition  7.1.  Let  F  be  an  ISOP  of  ST{n,l)-  In  the  graph 
representation  Gf  of  F 

1.  Gf  has  nodes  xi,X2, . . . ,  and  x„,  and 

2.  Gf  has  an  edge  from  Xi  to  Xj  iff  XiXj  is  a  PI  in  F. 
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TABLE  3 

Distribution  of  ISOPs  in  ST{n,  1)  Functions 


#  of  Pis 

5T(3,1) 

5T(4,1) 

5r(5, 1) 

5r(6,i) 

5T(7, 1) 

3 

2 

4 

3 

6 

5 

36 

24 

6 

16 

360 

120 

7 

560 

3600 

720 

8 

125 

13530 

37800 

9 

9270 

282660 

10 

1296 

435330 

11 

170352 

12 

16807 

Total 

5 

58 

1069 

27816 

943669 

Example  7.1.  Fig.  5  shows  the  graph  representations  of  the 
MSOP  and  WSOP  for  5^(3, 1)  (shown  in  Fig.  1). 

We  show  that  the  graph  representation  of  an  ISOP  of  F 
has  a  special  property. 

Definition  7.2.  A  directed  graph  G  is  strongly  connected  iff  for 
every  pair  of  vertices  (a,  b)  in  G,  there  is  a  path  from  atob  and 
from  b  to  a.  A  directed  graph  G  is  minimally  strongly 
connected  iff  it  is  strongly  connected  and  the  removal  of  any 
edge  causes  G  not  to  be  strongly  connected. 

Theorem  7.1.  Let  Gp  be  a  graph  representation  of  F.  F  is  an 
ISOP  of  ST{n,  1)  iff  Gp  is  minimally  strongly  connected. 

Proof.  See  the  Appendix.  □ 

The  graph  representations  of  the  MSOP  and  WSOP  of 
ST{3, 1),  shown  in  Fig.  5,  are  both  strongly  connected,  as 
they  should  be  by  Theorem  7.1.  Since  each  edge  represents  a 
prime  implicant,  an  MSOP  has  a  graph  representation  with 
the  fewest  edges.  This  observation  facilitates  the  enumera¬ 
tion  of  MSOPs. 


Theorem  7.2.  The  number  of  MSOPs  for  ST{n,  1)  is  {n  —  1)!. 
Proof.  See  the  Appendix.  □ 

The  graph  representation  allows  a  characterization  of 
ISOPs.  Specifically,  complementing  all  variables  in  an  ISOP 
of  ST{n,  1)  is  equivalent  to  reversing  the  direction  of  all 
edges  in  the  graph  representation  Gp  of  F.  If  Gp  is 
minimally  strongly  connected,  then  the  graph  obtained 
from  Gp  hy  reversing  the  direction  of  all  edges  is  also 
minimally  sfrongly  connected.  This  proves: 

Lemma  7.1.  If  F  is  an  ISOP  of  ST{n,  1),  then  the  SOP  derived 
from  F  by  complementing  all  variables  is  an  ISOP  of 

ST{n,l)- 

Example  7.2.  When  all  variables  are  complemented,  the 
graph  representations  of  the  ISOPs  shown  in  Fig.  5 
produce  the  graphs  in  Fig.  6,  which  also  represent  an 
MSOP  and  a  WSOP. 

It  is  important  to  note  the  difference  between  changing  an 
ISOP  F  and  changing  the  function  realized  by  F.  That  is,  an 
ST{n,  k)  function  is  unchanged  by  a  complementation  of 
variables,  i.e.  it  is  a  self-anti-dual  function  [32]. 

However,  an  ISOP  for  an  ST{n,  k)  function  may  or 
may  not  be  changed  when  all  variables  are  comple¬ 
mented.  For  example,  F  =  xix^k  X2X3  V  X3X1  is  an  ISOP 
for  S'T(3,  !)•  Complementing  all  variables  in  F  yields, 
F  =  X1X2  V  X2X3  V  x^xi,  a  different  ISOP. 

It  is  interesting  that  the  WSOP  for  ST{3, 1)  is  unchanged 
by  a  complementation  of  all  variables,  as  can  be  seen  by 
comparing  Fig.  5b  with  Fig.  6b.  The  invariance  of  an  ISOP 
with  respect  to  complementation  of  all  literals  is  a  unique 
characteristic  of  WSOPs,  as  shown  in  the  next  result. 

Lemma  7.2.  Let  F  be  an  ISOP  of  ST{n,  1).  F  is  a  WSOP  iff 
complementing  all  variables  in  F  leaves  F  unchanged. 

Proof.  See  the  Appendix.  □ 


TABLE  4 

Distribution  of  ISOPs  in  Arithmetic  Functions 


#  of  Pis 

INC6 

WGT5 

ROT6 

LOG5 

ADR3 

SQR5 

22 

13942125 

12 

23 

1088 

286 

24 

3072 

1010 

25 

2640 

1326 

26 

720 

478 

9452 

27 

64 

281252 

28 

3599288 

29 

20725014 

30 

58836676 

31 

59049 

447561 

92597382 

32 

64224 

83902808 

33 

449865 

42813004 

34 

64224 

11297310 

35 

2304 

1340364 

36 

54176 

37 

280 
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TABLE  5 

Number  of  Pis  Produced  by  Various  Algorithms  on  Three  Benchmark  Functions 


Algorithm 

SQR5 

SQR6 

5^(9, 3)** 

Quine-McCluskey  [22] 

26  MSOP 

47  MSOP 

84  MSOP 

ESPRESSO-IIC  [3] 

26  MSOP 

49 

85 

MINI-APL  [15] 

- 

49 

85 

POP-C  [9] 

- 

53 

148 

CAMP  [1] 

- 

- 

130 

Minato-Morreale  [19,  20] 

30 

58 

148 

WIRR* 

31 

71 

136 

Algorithm  5.1 

37  WSOP 

t 

t 

*  Algorithm  to  find  an  ISOP  with  many  Pis. 

**  5T(9.3)  is  equivalent  to  9SYM. 
f  Memory  overflow  precluded  generation  of  these  values. 

—  For  these  algorithms,  there  are  no  juiblished  results.  We 
do  not  have  the  code. 

Values  marked  MSOP  and  WSOP  are  known  to  be  MSOP 
and  WSOP. 


It  is  interesting  that  Lemma  7.2  does  not  generalize  to 
ST{n,  k).  Specifically,  for  S'T(5, 2),  the  ISOP 

F  =  X1X2XZX5  V  X3X5X1X2  V  3:10:30:2314  V  a;2a;4a:ia:3 

V  3;io;40:20:5  V  3;2a;5a:io:4  V  3;io;50:33:4  V  X3X4X1X5 

V  X2X3X4X3 V  X4X3X2X3 

is  invariant  with  respect  to  complementation  of  all 
variables.  However,  it  is  an  MSOP  and  not  a  WSOP. 

We  can  also  enumerate  WSOPs  as  follows: 

Theorem  7.3.  The  number  of  WSOPs  for  ST{n,  1)  is 
Proof.  See  the  Appendix.  □ 

The  graph  representation  allows  the  enumeration  of 
other  classes  of  ISOPs.  For  example,  we  can  enumerate 
ISOPs  that  have  one  more  PI  than  is  in  the  MSOP. 
Specifically, 

Theorem  7.4.  The  number  of  ISOPs  for  ST{n,  1)  with  n  +  1  Pis 


I  I  n  —  1 


Proof.  See  the  Appendix.  □ 

By  comparing  the  number  of  MSOPs  with  either  the 
number  of  WSOPs  or  the  number  of  ISOPs  with  one  more  PI 


than  in  the  MSOP,  we  find  that  the  former  is  much  less  than 
either  of  the  latter  for  large  n.  That  is,  as  n  approaches 
infinity,  the  ratio  of  MSOPs  to  WSOPs  approaches  0  (use 
Stirling's  formula  to  replace  (n  —  1)!  in  the  expression  for 
the  number  of  MSOPs).  This  proves  the  following: 
Theorem  7.5.  The  fraction  of  ISOPs  for  ST{n,  1)  that  are 

MSOPs  approaches  0  as  n  approaches  infinity. 

It  is  interesting  that  the  ratio  of  the  number  of  ISOPs  with 
n  +  1  Pis  (one  more  PI  than  is  in  an  MSOP)  to  the  number  of 
WSOPs  also  approaches  0  as  n  approaches  infinity.  This 
suggests  that  WSOPs  are  much  more  common  than  minimal 
or  near-minimal  ISOPs. 

8  Concluding  Remarks 

The  existence  of  an  algorithm  that  finds  the  worst  sum-of- 
products  expression  for  a  class  of  functions  is  surprising.  It 
counters  our  expectation  that  a  heuristic  algorithm  should 
perform  "reasonably"  well.  Also,  the  large  difference 
between  the  size  of  the  worst  and  the  best  expression  is 
especially  compelling  since  such  an  algorithm  will  perform 
very  poorly.  It  is,  therefore,  an  interesting  question  of 
whether  there  are  other  algorithms  and  other  functions  that 
exhibit  the  same  characteristics. 


XI  X] 


X3X2 


(a)  (b) 

Fig.  5.  Graph  representations  of  the  MSOP  and  WSOP  for  5T(3, 1). 
(a)  MSOP.  (b)  WSOP. 


X2X3 

(a)  (b) 

Fig.  6.  Graph  representations  of  Fig.  5  with  aii  variabies  compiemented. 
(a)  MSOP.  (b)  WSOP. 
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We  show  a  multiple-output  function  where  the  worst 
and  the  best  ISOPs  differ  greatly  in  size.  Specifically,  a 
decoder  with  2"  outputs  and  n  inputs  realizes  a  function 
where  a  WSOP  has  at  least  2"  Pis  and  an  MSOP  has  at  most 
2n  Pis.  Since  this  is  a  commonly  used  logic  function, 
disparity  in  the  size  of  WSOPs  and  MSOPs  carmot  be 
viewed  as  a  characteristic  of  contrived  functions  only. 

Although  computationally  intensive,  enumeration  of  the 
ISOPs  for  representative  functions  gives  needed  insight  into 
the  problem.  We  show  an  algorithm  to  compute  all  ISOPs  of 
a  given  function.  We  apply  it  to  benchmark  functions  and 
show  there  are  significant  differences  in  the  distributions  of 
ISOPs.  That  is,  some  functions  have  a  narrow  distribution, 
where  the  WSOP  is  nearly  or  exactly  the  same  size  as  the 
MSOP.  These  tend  to  be  easy  to  minimize.  For  example,  for 
unate  functions  [17]  and  parity  functions,  there  is  exactly 
one  ISOP.  Such  functions  are  classified  as  "trivial"  in  the 
Berkeley  PLA  Benchmark  Set  (e.g.,  ALUl,  BCD,  DIV3, 
CLPl,  C014,  MAX46,  NEWPLA2,  NEWBYTE,  NEWTAG, 
and  RYY6)  [26].  Other  functions  display  a  wide  range  and 
tend  to  be  hard  to  minimize.  For  example,  9SYM  or  SYM9 
(S'T(9,3))  has  a  wide  range,  i.e.,  the  number  of  Pis  in  a 
WSOP  and  an  MSOP  is  148  and  84  Pis,  respectively.  This 
function  is  known  to  be  hard  to  minimize. 

For  a  class  of  functions,  we  provide  an  analysis  showing 
that  the  number  of  MSOPs  is  significantly  smaller  than  the 
number  of  WSOPs.  That  is,  by  showing  a  correlation  with 
directed  graphs,  we  enumerate  all  MSOPs  and  all  WSOPs  of 
the  class  and  show  that  the  number  of  MSOPs  and  WSOPs 
is  (n  —  1)!  and  respectively.  As  n  increases,  the  ratio  of 
Pis  in  a  WSOP  to  the  Pis  in  an  MSOP  grows  without  bound. 
This  suggests  such  functions  are  hard  to  minimize. 

A  complete  understanding  of  the  minimization  process 
will  require  knowledge  of  the  search  space  and  how  various 
algorithms  progress  through  it.  However,  such  an  under¬ 
standing  is  not  likely  to  be  achieved  in  the  near  future.  Our 
research  suggests  that  there  is  merit  to  understanding  the 
correlation  between  the  degree  of  difficulty  in  minimizing  a 
function  and  the  distribution  of  its  ISOPs. 

Appendix 

Lemma  2.1. 


t{WSOP  :  5[oa,3.4.6,7})  =  70. 

Proof.  There  are  two  steps.  In  the  first  step,  we  prove  that 
an  ISOP  with  70  Pis  exists  for  this  function.  In  the  second 
step,  we  show  that  it  is  a  WSOP.  For  the  first  step,  it  is 
convenient  to  view  the  symmetric  function  as  having 
three  parts.  Specifically, 

‘^{0.1, 3, 4.6, 7}  =  ‘^{0,1}  ‘^{3,4}  ‘^{6,7}- 

A  WSOP  is  obtained  by  finding  a  WSOP  of  each  of  the 
three  parts  separately.  Consider  the  7-bit  Hamming  code 
shown  in  Table  6. 

For  each  code  word,  create  a  PI  that  covers  two 
minterms  by  replacing  one  of  the  most  abundant  bits  in 


TABLE  6 

A  7-Bit  Hamming  Code 
Code  word  Code  word 


0  0  0  0  0  0  0 
0  0  0  1  1  1  1 
0  0  1  0  0  1  1 
0  0  1110  0 
0  10  0  10  1 
0  10  10  10 
0  110  110 
0  1110  0  1 


1111111 
1  1  1  0  0  0  0 
110  110  0 
1  1  0  0  0  1  1 
10  110  10 
10  10  10  1 
1  0  0  1  0  0  1 
1  0  0  0  1  1  0 


the  code  word  by  a  don't  care.  In  the  case  of  code  word, 
0000000,  this  creates  seven  Pis,  each  of  which  covers  the 
minterm  with  all  variables  0  and  one  minterm  with 
exactly  one  1.  This  covers  all  minterms  of 
Similarly,  seven  Pis  generated  from  code  word  1111111 
cover  all  minterms  of  S'jg  ^ . 

All  of  the  remaining  14  code  words  have  either  four  Os 
and  three  Is  or  four  Is  and  three  Os.  For  each,  create  four 
Pis  by  changing  one  of  the  four  logic  values  in  the 
majority  to  a  don't  care.  Collectively,  the  four  Pis  cover 
the  original  code  word  and  four  words  that  are  a 
distance  one  away  from  the  code  word.  Because  the 
distance  between  any  pair  of  code  words  is  at  least  three, 
a  change  in  a  single  bit  of  a  code  word  in  the  Hamming 
code  creates  a  word  that  is  not  a  code  word  that  is 
distinct  from  either  another  code  word  or  a  word  that  is 
one  bit  different  from  another  code  word.  This  implies 
that  those  minterms  a  distance  one  away  from  a  code 
word  are  covered  by  at  most  one  PI.  It  follows  that  each 
PI  is  irredundant.  Since  each  of  the  14  code  words 
corresponds  to  a  set  of  Pis  that  cover  five  distinct 
minterms,  there  are  14  x  5  =  70  minterms  total.  On  the 
other  hand,  the  number  of  minterms  for  is 

(D  (D  =  35  -I-  35  =  70.  It  follows  that  these  Pis  cover 
all  the  minterms  of  5'|3 In  all,  7  -I-  56  -I-  7  =  70  Pis  cover 
all  of  the  8  -L  35  -L  35  -L  8  =  86  minterms  of  fhe  function. 
It  follows  that  this  set  of  Pis  is  a  cover  for  Sj^  13457}. 
Further,  it  is  an  irredundant  cover  and  we  have  an  ISOP. 

We  have  proven  that  an  ISOP  with  70  Pis  exists  for 
S'jo  1  3  4  6  7}-  show  that  this  is  a  WSOP  by  showing  that 
no  more  than  seven,  56,  and  seven  Pis  can  cover  the 
minterms  in  510  1},  <S'}3  4},  and  <S'{g7},  respectively.  Since 
S’jgi}  and  <S'|g7}  are  monotone  functions,  their  ISOPs  are 
unique.  Each  consists  of  seven  Pis.  For  <S'|3  4},  the  ISOP  of 
56  Pis  above  covers  the  70  minterms  associated  with  this 
function.  On  the  contrary,  assume  that  the  proposed 
ISOP  is  not  a  WSOP.  Thus,  there  is  a  set  of  p  >  56  Pis  that 
forms  an  ISOP  for  these  minterms.  Each  PI  covers  exactly 
two  minterms,  for  a  total  of  2p  >  112  instances  of  a  PI 
covering  a  minterm.  Let  m\  and  m>i  be  the  number  of 
minterms  covered  by  one  and  more  than  one  PI, 


944 


IEEE  TRANSACTIONS  ON  COMPUTERS,  VOL.  50,  NO.  9,  SEPTEMBER  2001 


respectively.  m>i  =  70  — mi.  Since  the  set  of  Pis  is 
irredundant,  each  PI  covers  at  least  one  minterm  that  is 
not  covered  by  any  other  PI.  Thus,  mi  >  p  >  56.  It 
follows  that  2p  +  3mi  >  280.  Further,  2p  —  mi  >  280  — 
4mi  =  4(70  —  mi)  and  we  can  write 


2p  —  mi 

70  —  mi 


>  4. 


(A.l) 


Here,  the  numerator  is  the  number  of  instances  in  which 
a  PI  covers  a  minterm  that  is  covered  by  more  than  one 
PI,  while  the  denominator  represents  the  number  of 
minterms  covered  by  more  than  one  PL  Since  this  ratio 
exceeds  four,  by  the  Pigeonhole  Principle,  there  is  at  least 
one  minterm  covered  by  at  least  five  Pis.  But,  this  is 
impossible;  each  minterm  is  covered  by  no  more  than 
four  Pis  (i.e.,  each  code  word  is  covered  by  a  PI  derived 
from  a  code  word  by  converting  one  of  the  four  most 
abundant  variables,  0  or  1,  to  a  don't  care).  Thus,  it  must 
be  that  the  proposed  ISOP  is  a  WSOP.  □ 


Theorem  3.1. 

2.  T{MSOP:ST{n,k))={l). 

3.  t{WSOP  :  ST{n,  k))  >  2(^)  -  f*). 

Proof. 

1.  An  implicant  of  ST{n,  k)  has  the  form 


^  •  •  •  •  ^  •  'T* .  .  •  « 


That  is,  for  this  implicant  to  be  1,  at  least 
k  variables  {xi^ ,  ,  •  •  • ,  and  XiJ  must  be  1  and  at 

least  k  variables  must  be  0  ,  •  •  • ,  and 

XiJ,  where  2k  <  n.  This  implicant  is  prime; 
deleting  a  literal  creates  an  implicant  that  is  1 
when  ST{n,  k)  should  be  0.  Specifically,  deleting 
Xi^,Xi^,. . . ,  or  Xii^  creates  an  implicant  that  is  1 
when  less  than  n  variables  are  1,  while  deleting 
^i„-k+i  I  ^in-k+2 1  •  •  •  /  or  creates  an  implicant  that  is 
1  when  more  than  n  —  k  variables  are  1. 

The  number  of  such  Pis  is  the  number  of  ways 
to  separate  n  variables  into  three  parts,  where 
order  within  a  part  is  not  important.  This  is  the 


multinomial  =  k\(n%.k\  ■ 

First,  we  show  that  (^)  is  a  lower  bound  on  the 
number  of  Pis  of  ST{n,  k).  Then,  we  show  a  set  If 
of  (^‘)  Pis  that  covers  all  and  only  minterms  in  the 
function  ST{n,  k).  It  follows  that  the  OR  of  all  Pis 
in  n  is  an  MSOP  for  ST(n,  k). 

Consider  the  set  MI  of  (^)  minterms  of 
the  form  Xi^Xi^  ■  ■  ■  where  ij  e 

{1, 2, . . . ,  n  —  1}  and  n  >  2k.  That  is,  MI  consists 
of  minterms  that  are  1  when  exactly  k  of  the  n 
variables  are  0.  No  PI  for  ST{n,  k)  covers  two  or 
more  minterms  in  MI.  As  such,  MI  is  a  set  of 
independent  minterms  and  at  least  (^)  Pis  are 
needed. 


n  is  formed  as  follows:  For  each  minterm  mt  in 
MI,  apply  Algorithm  1.1  below,  producing  P,nt,  a 
PI  that  covers  mt.  Add  Pmt  to  TI.  Since  no  PI 
covers  two  or  more  minterms  in  MI,  n  has  (^) 
distinct  Pis.  Since  Pmt  has  exactly  k  Os  and  fc  Is,  it 
covers  only  minterms  in  ST{n,  k).  Next,  we  show 
that  n  covers  all  minterms  in  ST{n,  k)  by  applying 
Algorithm  1.1  to  an  arbitrary  minterm  mt'  of 
ST{n,  fc),  producing  Pmt',  a  PI  that  covers  mt'. 
Pmf  e  n,  as  follows:  Form  a  minterm  mt"  in  MI 
from  Pmt’  by  setting  all  -s  to  Is.  Applying 
Algorithm  1.1  to  mt"  yields  Pmt"  that  is  identical 
to  Pmt',  from  which  we  can  conclude  Pmt'  £  If. 

Algorithm  1.1  (Produce  a  PI  that  covers  a  given 
minterm). 

Input:  Minterm  mt  =  mt{0)mt{l) . . .  mt{n  —  1) 
Output:  Prime  implicant 

Pmt  =  P’mt(0)Pmt(l)  .  .  .  Pmt{n  -  1) 

(Initially,  Pmtii)  =  “  for  all  i  where  0  <  f  <  n  —  1) 

1 .  ZeroOnePairs  =  0 

2.  Repeat  until  ZeroOnePairs  =  fc  do  j 
if 

{Pmt{i)Pmt{i  +  1)  •  •  •  Pmt{i  +  s)  £  1}  “ 

and  mt{i)mt{i  +  s)  =  01) 

then 

{_Pmt(i)Pmt{i  +  s)  <—  01 

and  ZeroOnePairs  <—  ZeroOnePairs  +1}, 

where  index  addition  is  mod  n  (in  this 
algorithm,  we  assume  that  subscript  indices 
range  from  0  to  n  —  1,  i.e.,  the  variables  are 
xo,Xi,  -  ■  ■ ,  x„-i)  and  where  1  <  s  <  n  —  1. 

It  is  straightforward  to  show  that  Algorithm  1.1 
produces  a  PI  in  11  if  the  minterm  input  has  at 
least  fc  Is  and  fc  Os.  This  PI  may  be  different 
depending  on  the  values  of  i  chosen  in  each 
repetitive  step. 

3.  Assume  t(WSOP  :  ST(n  -  1,  fc))  >  2("-i)  -  (f). 
We  can  form  an  ISOP  of  ST{n,  fc)  as  follows: 

F  =  x„Fi  V  x„F2  V  F3, 

where 

a.  Fi  is  an  SOP  such  that  each  product  term 
is  formed  as  the  AND  of  i)  a  set  Xi  of 
fc— 1  complemented  variables,  where 
Xi  C  X  —  {x„},  and  of  ii)  fc  uncomplemented 
variables  from  X  —  {xn}  —  Xi,  where  the 
indices  of  fhe  uncomplemented  variables 
are  all  as  small  as  possible  (given  the  choice 
of  Ai).  Because  Xi  can  be  chosen  in 
ways,  Fi  has  product  terms. 

b.  F2  is  an  SOP  such  that  each  product 
term  is  formed  as  the  AND  of  i)  a  set 
X2  of  fc  —  1  uncomplemented  variables. 
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where  X2  C  X  —  {xn}  and  of  ii)  k  comple¬ 
mented  variables  from  X  —  {x„}  —  X2,  where 
the  indices  of  the  complemented  variables  are 
all  as  small  as  possible.  Because  X2  can  be 
chosen  in  Qr])  ways,  F2  has  product 

terms. 

c.  F3  is  one  of  the  WSOPs  for  ST{n—l,k). 
From  the  inductive  hypothesis,  F3  has 

-  Ck)  Pis- 

F  is  an  expression  for  ST{n,  k)  as  follows: 
Consider  a  minterm  m  in  ST{n,  k).  If  m  has  at 
least  k  Is  and  at  least  k  Os,  regardless  of  the 
value  of  x„,  then  m  is  covered  by  F^.  It  m  has 
exactly  k  Os,  including  a  0  value  for  Xn,  and  at 
least  k  Is,  then  it  is  covered  by  x„Fi.  If  m  has 
exactly  k  Is,  including  a  1  value  for  a;„,  and  at 
least  k  Os,  it  is  covered  by  XnF2.  Thus,  F 
covers  all  minterms  in  ST{n,  k).  Because  each 
PI  in  XnF\  and  XnF2  (and  also  in  F3)  has 
exactly  k  uncomplemented  and  exactly 
k  complemented  variables,  F  covers  only 
minterms  with  at  least  k  Os  and  at  least  A:  Is, 
i.e.,  only  minterms  in  ST{n,  k).  It  follows  that 
F  is  an  SOP  for  ST{n,k). 

If  F  has  no  redundant  Pis,  its  (^“J)  -I- 
iT-l)  +  -  (?)  or  2(3  -  (^3  Pis  form 

an  ISOP  for  ST{n,  k).  Thus, 

r{WSOP  :  ST{n,k))  >  2Q  - 

Next,  we  show  that  F  has  no  redundant 
Pis.  First,  each  PI  in  XnFi  covers  a  minterm  m! 
having  k  Os  and  n  —  fc  Is  that  is  not  covered 
by  any  combination  of  Pis  from  XnF2  and  F^. 
Thus,  no  PI  in  XnFi  is  redundant.  By  a  similar 
argument,  no  PI  in  XnF2  is  redundant. 
Second,  no  PI  in  F3  is  redundant,  as  follows: 
Since  F3  is  a  WSOP  for  ST{n  —  1,  k),  no  PI  in 
F3  is  covered  by  the  OR  of  one  or  more  Pis  in 
F3.  If  the  OR  of  Pis  from  XnFi  and  x„F2  covers 
a  PI  P  from  P3,  then  it  follows  that  at  least 
one  product  term  Pi  in  Pi,  when  ANDed 
with  at  least  one  product  term  P2  in  P2,  yields 
a  non-0  result.  Let 

Pi  —  Xsi  Xs2  .  •  ■  ^t2  •  •  ■ 

and 

P2  ^U2  •  •  ■  T-Ui  ^V2  •  *  •  ■ 

If  P1P2  ^  0,  no  Si  is  the  same  as  a  Uj  and  no  tp 
is  the  same  as  a  w,.  But,  ti,t2, ,  and  tk  were 
chosen  to  be  as  small  as  possible  without 
overlapping  si,s2,...,  and  st-i,  while 
Vi,V2,...,  and  Vk  were  chosen  to  be  as 
small  as  possible  without  overlapping 
Ui,U2,...,  and  Uk-i-  Consider  the  indices 
/={l,2,...,fc}.  The  smallest  index  in  I  that 
appears  neither  in  S'  =  {si,  S2,  •  ■  • ,  Sfc-i}  nor  in 
U  =  {ui,U2, . . . ,  Mjt-i}  appears  in  both  T  = 


{ti,t2,---,4}  and  V  =  {vi,V2, .  ■  ■  ,Vk},  caus¬ 
ing  P1P2  =  0,  a  contradiction.  □ 

Theorem  3.5.  ST{m,  ky  has  the  following  properties: 

1.  riCSOP  :  ST(m,  kf)  = 

2.  t{MSOP  :  ST{m,  kf)  =  (7)’'. 

3.  ywsop  :  ST(m,  kf)  >  [2(™)  -  y^)Y. 

Proof.  Items  1,  2,  and  3  follow  from  the  observation  that  a 
minterm  in  the  product  function  ST{m,ky  can  be 
viewed  as  the  AND  of  a  minterm  from  each  of  the  factor 
functions  ST{m,k).  Thus,  a  PI  of  ST{m,ky  can  be 
viewed  as  the  product  of  a  PI  from  each  ST{m,  k)  and 
Item  1  follows  directly. 

Also, 


t{MSOP  :  ST{m,ky)  < 

follows  directly.  That  is,  (^Y  is  an  upper  bound  on  the 
number  of  Pis  in  an  MSOP  for  ST{m,  ky,  as  an  ISOP  for 
ST{m,ky  can  be  formed  as  the  AND  of  Pis  from  the 
MSOPs  of  ST{m,  k) .  As  is  shown  by  Voight  and  Wegener 
[36],  certain  product  functions  can  have  fewer  Pis  in  their 
MSOPs  than  the  product  of  the  number  of  Pis  in  the 
MSOPs  of  the  factor  functions.  However,  when  the  factor 
functions  are  ST{m,  k),  we  can  observe  the  following:  Let 
M  be  the  set  of  minterms  covered  by  ST{m,  k),  in  which 
exactly  k  variables  are  uncomplemented.  Since  the  Pis  of 
ST{m,k)  have  exactly  k  uncomplemented  and  k  com¬ 
plemented  variables,  none  cover  two  or  more  minterms 
in  M.  It  follows  that  at  least  (™)  Pis  are  needed  to  cover 
ST{m,k).  It  follows  that  at  least  (ffY  Pis  are  needed  to 
cover  the  minterms  in  ST{m,  ky  that  are  the  product  of 
minterms  in  the  set  M  for  each  factor  function.  Thus, 

t{MSOP  :  ST{m,ky)  =  . 

A  WSOP  for  ST{m,  ky  can  be  formed  as  the  product  of 
WSOPs  for  each  factor  function.  Since  2(™)  —  is  a 
lower  bound  on  the  number  of  Pis  in  each  factor 
function.  Item  3  follows  directly.  □ 

Theorem  7.1.  Let  Gp  be  the  graph  representation  of  F.  F  is  an 
ISOP  of  ST {n,  1)  iff  Gp  is  minimally  strongly  connected. 

Proof,  (if)  Consider  a  minimally  strongly  connected  digraph 
Gp,  where  F  is  its  corresponding  SOP.  Thus,  for  every 
edge  {xj,  Xi)  in  Gp,  there  is  an  implicant  XjXi  in  F.  On  the 
contrary,  assume  that  F  is  not  an  ISOP  for  ST{n,  1).  That 
is,  either  1)  F  does  not  cover  ST{n,  1)  or  2)  F  covers 
ST{n,  1)  but  has  a  redundant  PL 

Consider  1).  If  F  does  not  cover  ST{n,  1),  then  there  is 
a  minterm  mt  such  that  mt  either  a)  has  no  complemen¬ 
ted  variables  or  no  uncomplemented  variables,  but  F  is  1 
for  this  assignment  or  b)  has  at  least  one  uncomplemen¬ 
ted  variable  and  at  least  one  uncomplemented  variable, 
but  F  is  0  for  this  assignment.  The  first  part,  a),  is  not 
possible;  all  Pis  cover  only  minterms  with  at  least  one 
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complemented  variable  and  at  least  one  uncomplemen¬ 
ted  variable.  The  second  part,  b),  is  also  not  possible,  as 
follows:  Because  Gp  is  strongly  connected,  there  is  a  path 
Xj  =  Xki,Xk2,  -  ■  ■ ,  Xk^  =  Xi,  from  Xj  to  Xi.  Since  Xj  =  0  and 
Xi  =  1,  the  assignment  of  values  to  Xk^Xk^,  ■  ■  ■ ,  and  Xk„ 
corresponding  to  mt  has  the  property  that  there  is  an  s 
such  that  Xk^  =  0  and  Xk,+i  =  1-  The  corresponding  PI 
Xk,Xk,+i  is  in  F  and  is  1  for  the  assignment  associated 
with  mt.  Thus,  F  covers  mt. 

Consider  2).  If  F  covers  ST{n,  1),  but  has  a  redundant 
PI,  XjXi,  then  F',  which  is  F  with  xjXi  removed,  also 
covers  ST{n,l).  But,  Gp'  is  Gp  with  one  edge,  (xjXi), 
removed.  It  follows  that  Gp  is  not  minimally  strongly 
connected. 

(only  if)  Let  Gphea  graph  representation  of  an  ISOP  F 
of  S'T(n,  1).  Assume,  on  the  contrary,  that  Gp  is  not 
minimally  strongly  connected.  That  is,  either  1)  Gp  is  not 
strongly  cormected  or  2)  Gp  is  strongly  connected,  but  is 
not  minimal. 

If  Gp  is  not  strongly  cormected,  there  are  two 
nodes  Xj  and  Xi  such  that  no  path  exists  from  Xj  to 
Xi.  Let  Suc[xj)  be  the  set  of  all  nodes  for  which 
there  is  a  path  from  Xj,  i.e.,  all  successors  of  Xj.  Let 
Pre{xi)  be  the  set  of  all  nodes  for  which  there  exists 
a  path  to  Xi,  i.e.,  all  predecessors  of  a;,.  Consider  a 
minterm  mt  that  is  0  for  Xj  and  all  variables  associated 
with  nodes  in  Suc{xj)  and  is  1  for  Xi  and  all  variables 
associated  with  nodes  in  Pre{xi).  Since  there  is  no  path 
from  Xj  to  Xi,  Suc{xj)  n  Pre{xi)  =  f  and  such  an 
assignment  assigns  exactly  one  value  to  nodes  in 
Suc{xj)yj  Pre{xi)\J  {xj^Xi}.  Choose  the  values  of  all 
other  variables  to  be  1  (or  0).  No  edge  in  Gf  has  a  0  at 
its  tail  and  a  1  at  its  head.  Thus,  all  Pis  are  0  and  F  is  not 
an  SOP  for  ST{n,  1).  It  is,  thus,  not  an  ISOP,  a  contra¬ 
diction. 

If  Gp  is  strongly  cormected,  but  is  not  minimal,  there 
is  at  least  one  edge  {xj,  Xi)  that  can  be  removed  without 
affecting  the  connectedness  of  Gp.  It  follows  that  Gpi, 
where  F'  is  F  with  XjXi  removed,  is  a  graph  representa¬ 
tion  of  F' ,  an  SOP  for  the  same  function  as  F.  Thus,  F  is 
an  SOP,  but  not  an  ISOP,  contradicting  the  assumption.D 

Theorem  7.2.  The  number  of  MSOPs  for  ST{n,  1)  is  {n  —  1)!. 

Proof.  From  Theorem  7.1,  an  MSOP  for  ST{n,  1)  corre¬ 
sponds  to  a  directed  graph  with  the  fewest  edges  which 
is  strongly  connected,  but  is  not  strongly  connected 
when  any  edge  is  removed.  Such  a  graph  is  a  directed 
cycle  of  arcs  through  all  variables.  As  such,  it  represents 
a  cyclic  permutation  on  the  variables.  The  number  of 
such  permutations  is  (n  —  1)!.  □ 

Lemma  7.2.  Let  F  be  an  ISOP  of  ST{n,l)-  F  is  WSOP  iff 
complementing  all  variables  in  F  leaves  F  unchanged. 

Proof,  (if)  Let  F  be  an  ISOP  that  is  unchanged  by  a 
complementation  of  all  variables.  This  implies  that  if  XiXj 
is  a  PI  of  F,  then  so  also  is  XjXi.  It  follows  that  if  {xi,Xj)  is 
an  edge  in  Gp,  then  {xj,Xi)  is  an  edge  in  Gp.  That  is,  all 
edges  between  nodes  occur  in  pairs,  one  going  one  way 
and  the  other  going  the  other  way.  In  such  a  graph,  there 
are  n  —  1  pairs  or  2n  —  2  edges  in  all.  (Replace  each  pair 
by  an  undirected  edge.  If  the  directed  graph  is  strongly 


connected,  the  undirected  graph  must  be  connected. 
From  Harary  [14],  there  must  be  n  —  1  edges.)  Since  there 
are  2n  —  2  edges  in  Gp,  there  are  2n  —  2  Pis  in  F  and, 
thus,  F  is  a  WSOP. 

(only  if)  Let  F  be  a  WSOP  of  ST(n,  1).  We  show  that 
Gp  consists  of  cycles  of  length  2  only.  Thus,  if  XiXj  is  a  PI 
of  F,  so  also  is  XjXi.  It  follows  that  complementing  all 
variables  of  F  leaves  F  unchanged.  Suppose  that  Gp 
contains  a  cycle  of  length  m,  where  m  >  2.  Such  a  cycle 
represents  a  strongly  connected  subgraph  of  Gf  in  which 
there  are  m  edges.  However,  the  cycle  can  be  replaced  by 
a  minimally  strongly  connected  graph  with  more  edges 
(e.g.,  where  all  edges  occur  in  pairs).  The  result  is  a 
strongly  cormected  graph,  where  the  deletion  of  an  edge 
leaves  it  unconnected,  which  has  more  edges  than  the 
original  version.  This  contradicts  the  statement  that  F  is 
a  WSOP.  □ 

Theorem  7.3.  The  number  of  WSOPsfor  ST{n,  1)  is  rP~^. 

Proof.  From  Theorem  7.1,  a  WSOP  for  ST{n,  1)  corresponds 
to  a  minimally  strongly  connected  graph  with  the  largest 
number  of  edges.  We  show  that  this  graph  consists  of 
cycles  of  length  2  exclusively,  as  follows:  Suppose,  on  the 
contrary,  the  graph  has  a  cycle  of  length  m  >  2.  There  are 
m  edges  in  this  cycle.  However,  this  subgraph  can  be 
replaced  by  a  subgraph  with  more  edges,  2(m  —  1).  It 
follows  that  the  original  graph  does  not  represent  a 
WSOP. 

Each  cycle  of  length  2  cormects  two  nodes  by  edges  in 
the  two  directions.  Replace  each  pair  of  edges  by  an 
undirected  edge,  forming  an  undirected  tree  with  n  —  1 
edges.  Thus,  there  are  2(n  —  1)  Pis  in  a  WSOP  for 
ST{n,l)-  This  proves  Theorem  3.2.  It  follows  that  the 
number  of  WSOPs  is  the  number  of  undirected  trees  on  n 
labeled  nodes.  Cayley  [5]  in  1889  showed  that  this 
number  is  rP~^.  □ 

Theorem  7.4.  The  number  of  ISOPs  for  ST{n,  1)  with  n  -|- 1  Pis 
is 


Proof.  From  Theorem  7.1,  an  ISOP  of  ST{n,  1)  that  has  n  -I- 1 
Pis  corresponds  to  a  minimally  strongly  connected  graph 
with  n  -I- 1  edges.  All  graphs  with  this  property  have  two 
cycles  of  nodes,  of  which  i  are  common,  where 
1  <  i  <  n  —  2.  We  can  represent  each  instance  as  a 
permutation  of  nodes  that  has  been  divided  into  three 
nonempty  sets,  the  i  common  nodes,  nodes  A^i,  in  one 
cycle  only,  and  nodes  N2  in  the  other  cycle  only.  Ni  and 
N2  must  be  nonempty  since  an  empty  set  corresponds  to 
a  redundant  edge.  There  are  nl  ways  to  permute  the 
n  nodes  and  ('‘“^)  ways  to  divide  them  into  three 
nonempty  sets.  However,  this  double  counts  graphs 
since  interchanging  Ni,  and  N2  does  not  change  the 
graph.  Thus,  the  total  number  of  graphs  with  n  +  I  edges 

is  !("->!•  □ 
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